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ABSTRACT 



This report lays the groundwork for further development 
of a microcomputer that provides full celestial navigation 
capability. The physical design specifications and design 
philosophy are investigated. The fix algorithm which 
allows determination of position without the use of H.O. 
Publications is implemented in the BASIC language. 

The application of microcomputer technology to devices 
of this nature brings to light the need for an information 
distribution system. Thoughts on a system of this purpose 
are presented. 
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I. INTRODUCTION 



This paper deals with two interleaved subjects. The 
two subjects will be discussed with varying levels of 
detail and scope. The more specific subject of the two is 
the conceptual design of a celestial navigation computer. 

A general approach toward a digital information distribu- 
tion system is the second subject. The celestial navigation 
computer is discussed first, as it is an element of the 
entire system, and the design specifics of this device 
relate better to the entire system than vice versa. 

The celestial navigation computer provides for manual 
input of information as well as documentation of the out- 
put for legal records. The main reason why a computer is 
beneficial for celestial positioning is that it eliminates 
errors made in recording, computing and plotting of celes- 
tial information. The design efforts attempt to eliminate 
manual recording of substantive information and manual 
plotting of lines of position (LOP's) for use in determining 
position. Automatic information recording is provided. 

This paper provides an analysis of the programming 
language, program I/O requirements and physiological require- 
ments of the device. The initial purpose of this paper was 
to design the device, then to see how this philosophy and 
design process could be applied to other building block 
devices applicable to the information distribution system. 
Comments to this point will be made in Section VI. 
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II. NEED FOR A CELESTIAL MICROCOMPUTER PERIPHERAL 



It would seem that with all of the fancy new gadgetry 
on the market today, the age-old method of navigating by 
the stars might just start to fade away. It is conceivable 
that this could happen in the civilian world; however, the 
military is often more demanding in the constraints placed 
upon its ships. The electronic devices which use electro- 
magnetic radiation in the position determination loop are 
subject to deception, jamming and detection (as in the 
case of radar), and support failure, e.g., destruction of 
navigational satellites or LORAN stations, etc. The only 
other navigation device not subject to standoff disabling 
is the ship's inertial navigation system. It must be 
pointed out that an inertial navigation system is only as 
good as its platform's stability and the stability only 
as good as its positional estimates. Thus the longer the 
platform goes without being updated, the worse the platform 
stability becomes- -truly a spiraling decline in accuracy. 

With all of the electronic fixing devices out of the 
picture, visual sighting or sonar positioning is all that 
is left. Visual sighting is passive and includes celestial 
position determination, and this is the reason for main- 
taining excellence in position determination by the stars. 

The microcomputer is a newcomer to this field. This 
technology properly applied allows for increased speed and 
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accuracy and the reduction of human errors. A need to 
determine position faster simply goes along with the increase 
in ship speed, e.g., new hydrofoil ships can conceivably 
reach speeds in excess of 60 kiaots. The need for increased 
accuracy is an economic as well as operational benefit. 
Accurate positioning is needed for proper and efficient 
deployment of weapons, efficient navigation and maneuvering. 

The submarine force has implemented a celestial naviga- 
tion system in the most complete terms. Completely computer 
controlled, several "spot shots" are obtained in a moment 
and the submarine's navigation systems can then be updated. 
The Air Force has used the ASTRO-TRACKER to aid in celestial 
navigation . 

In January, 1975, Carl Nuese and Dr. Alan Schneider 
provided the algorithms necessary to create a celestial 
navigation calculator using the Hewlett Packard Model 65 
programmable microcomputer [Ref. 8] . The peripheral device 
under design in this paper was intended to hit the mark 
somewhere between the HP-65 calculator and the submarine 
celestial navigation computer- -much closer to the HP-65 
side of the spectrum. 

The HP-65 calculator solution to navigation problem 
lacks the same sort of smoothness that placing a square 
block into a round hole does; with a lot of effort it works 
but it leaves a lot to be desired. The HP-65 general purpose 
calculator uses programmable magnetic cards to hold prestored 
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programs that are executable on the machine. Though the 
device is relatively inexpensive the following drawbacks 
are noted: 

1. Storage and maintenance of the program cards are 
cumbersome . 

2. I/O is nonexistent. 

3. Data format is not "conventional.” 

4. Record keeping must still be done manually. 

5. The device is subject to pilferage. 

6. Expansion of programs is difficult in its 100- 
instruction memory. 

The advantages are the following: 

1. Self contained (excluding the magnetic cards). 

2. Easier to operate than most terminals. 

3. Inexpensive. 

4. Good supporting vendor. 

5. Applicable to many other mathematical problems. 
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III. DESIGN PHILOSOPHY 



A. OBJECTIVES 

The objectives of the design of the device are simple: 

1. Implement Nuese and Schneider's algorithms in a 
device that does not have the disadvantages of the HP-65. 

2. Maintain the advantages of the HP-65 without appre- 
ciably increasing the cost/performance factor. 

3. Provide for integration into the overall information 
distribution system. 

4. Provide the fail safe/soft features and hardcopy 
printout capabilities. 

Two definitions will be stated at this time. 

Microcomputer: A small general purpose computer. It 

contains a microprocessor and a small number of (sup- 
port) LSI circuits. Usually mounted on one or two 
printed circuit boards. May or may not be programmable. 

Minicomputer: A computer with larger word sizes and 

faster execution times than a microcomputer. 

The proposal to specify a microcomputer for the celestial 
computer was based on the further observations that the micro- 
computer can be exactly specified as to the I/O structure, 
memory structure and contents, po\^fer supply designs, etc., 
whereas the minicomputer requires adjusting of all of the 
above mentioned properties to "fit" the minicomputer. The 
cost for a microcomputer is also much less than that of the 
minicomputer . 

An eight-bit-wide microcomputer consisting of a micro- 
processor, a 128 X 8 bit RAM, a 1024 x 8 ROM and peripheral 
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adapter cost around $473 in 1974, and by 1976 is estimated 
to cost under $100 [Ref. 6] . Though larger than the HP-65 
it is still less than ten times as large. An examination 
of data equipment suppliers will show that a small mini- 
computer will cost around $5,000 and higher. 

The device must be self-contained (except for power) , 
easy to operate and relatively inexpensive. The I/O struc- 
ture must allow for keyboard entry, and the formats for data 
must be naturalized (36° 15' 33" vice 36.1533°) to the prob- 
lem. Note that the algorithm expressed in Appendix A uses 
the latter form due to limitations of the BASIC language 
I/O structure. As the device would be a small production 
lot, vendor support would be minimal or nonexistent; yet 
the simplicity of the device should allow for easy onboard 
maintenance . 

B. PHYSICAL DESIGN SPECIFICS 

The device must be housed in one unit. The unit must 
allow for connection to power and to the data distribution 
system. The unit will have a simple ten-digit keyboard with 
a N/E (north or east) and S/W (south or west) key. A 16 
element five-by-seven dot LED must be placed above the key- 
board for verification of entries as the line printer prints 
a completed entry. 

A 30 character wide line printer (for printing instruc- 
tion, entries and calculation results) must be placed in 
the unit. The unit is illustrated in Figure 1, and a 
sample output strip is displayed in Appendix B. 
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Figure l~Celestial Navigation Microcomputer 
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IV. THE CELESTIAL FIX ALGORITHM 



The purpose o£ the proposed navigational computer is 
to enhance the speed, accuracy and ease with which calcula- 
tions such as those used in celestial navigation can be made 
The program presented in Appendix A removes the need for the 
user to depend on the Nautical Almanac [Ref. 13] and other 
required H.O. publications by computing the values supplied 
by these tables. 

When taking a fix several phases of sight reduction are 
considered: insertion of user values, correction of the 

sextant altitude, establishment of the position of the 
observed celestial body, solution of the navigational 
triangle and displaying the results. 

User values include Greenwich date, day, month and year, 
the observer's assumed position, the name of the body ob- 
served, the time (Greenwich Mean Time), the sextant altitude 
the observer's height of eye or dip and the sextant's index 
of correction. 

Correction of the sextant's altitude includes changes 
in altitude due to refraction, semi -diameter of the sun and 
moon and parallax for the moon. 

Establishment of the position of the observed celestial 
body on the earth (GP) involves the calculation of the 
Greenwich Hour Angle (GHA) and declination (Dec) of the body 
(see Figure 2) . Associated with each star are eight 



14 















■; A-> 






I 




♦ 





Figure 2-Geographic Position (GP) oi a Star 
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intermediate constants which are calculated and combined 
with universal constants (rate of precession of the equinox, 
rotation of the earth on its axis and aberration) to produce 
the GHA and declination. 

Constants which apply to all stars are the following: 

1. Rate of rotation of the earth with respect to the 
vernal equinox [Ref. 11] . 

2. The longitude rate of the sun based on the time 
between successive crossings by the sun of the vernal 
equinox [Ref. 10] . 

3. The obliquity for epoch 1975 [Ref. 10]. 

4. The annual precession rates for epoch 1975 in right 
ascension and declination [Ref. 10] . 

With these common constants, particular intermediate 
constants associated with each star are developed. Utiliz- 
ing the apparent right ascension and declination as found 
in Apparent Places of the Fundamental Stars [Ref. 14] , 
right ascension and declination can be computed to a 
reference time through the effects of precession and 
aberration. Other constants include the rate of increase 
of the mean sun longitude in degrees per day, precession 
in right ascension and declination in degrees per day, true 
right ascension and declination in degrees, and the time 
behavior of the corrections for aberration [Ref. 8] . 

When the intermediate constants are developed and 
applied to the algorithms, the computed GHA and declination 
of a star will result. Discussion of the computation of 
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the constants, GHA and declination is presented in the 
next section. With the algorithm for computing GHA and 
declination of a star implemented it now becomes a matter 
of selecting those stars necessary to perform a sight 
reduction yielding a navigational fix. 

Considering a three-star fix, the GHA and declination 
for each star is obtained as discussed above. Working with 
pairs of stars the solution of the associated navigational 
triangles will produce intersection points on the circles 
of equal altitudes (see Figure 3) . 

Three triangles are used to calculate a possible position 
as shown in Figure 4. The first triangle has vertices N 
(north pole), A and B. The angle A is either G1-G2 or 
360° - (G1-G2) . This angle represents the difference in the 
GHA's. Given the angle A and cosides dl and d2 (the 
declinations) , the law of cosines is applied to solve the 
third side t. Given the three cosides the angle B opposite 
coside dl is calculated. Given three cosides in the second 
triangle, HI, H2 and t, the angle B3 opposite HI is found. 

Letting m be the sign of the sin(Gl-G2), the angles 
B(l)=mB-B3 and B(2)=mB+B3 are found. Given angle B(l) and 
cosides d2 and H2, the third coside L(l) is found. Given 
the three cosides d2, H2 and L(l) the angle A(l) opposite 
H2 is found. In the same manner L(2) and A(2) are found. 

L(i) is equivalent to the latitude and A(i) is equivalent 
to the longitude of an intersection of the circles of equal 
altitude . 
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Figure 3-Navigational Fix from Altitudes 
of Tv?o Celestial Bodies 
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Fig^ire Triangles used for Navigational Fix 
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Repeating the algorithm considering the third celestial 
body and one of the previously used bodies, a second pair 
of positions will be found. A third repetition yields two 
more positions totaling six (the maximum) possible positions, 
as shown in Figure 5. 

These six positions are then investigated to determine 
the closest three by determining the shortest distance 
connecting these three points. With this resulting triangle 
the location of the two endpoints of a line which is one- 
third away from any side of the triangle is computed. The 
center of this line containing the two endpoints is the 
center of mass of the triangle and therefore is considered 
the celestial fix. 

The program will provide for accuracy in mathematical 
manipulation of numbers to one part in four billion using 
the BASIC interpreter at the Naval Postgraduate School and 
implementing the program on the INTEL 8080. This will pro- 
vide for the accuracies of .1 or .2 miles as presented in 
Ref. 8. 
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V. PROGRAM TMEORY AND IMPLEMENTATION 



A. THEORY 

When computing the constants needed for the program 
a common reference time is necessary to provide for the 
accuracy of results. If the data is referenced to the mean 
equinox and equator of date, or to a fixed epoch then 
no error is introduced. The reference time employed in 
the algorithms is 0000 29 Nov 1971 as suggested by Ref. 8. 
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where t is the present time measured in days and decimal 
parts of a day from the time t^ , a time in the recent past. 
Because it will later prove efficient to refer all events 
to an earlier time (t^ = 0000 29 Nov 1971) this leads to 

GHA = a(t - t^) - a(t - t^) + - a(t) (5) 

= a(t - t^) H- GHA^^ - a(t) (6) 

where GHA = GHA , - a(t, - t ) (7) 

yo yl 1 

The apparent right ascension a at time t is related to 
its known true value a' at some previous time through 
the effect of precession and aberration. 

Therefore (see "Supplement," [Ref. 11, pp , 39-49]) 

a(t) = a'(t^) + a(t - tj^) + Aa(t) (8) 

where 

a'(tj) = the true right ascension at time t^ 

a = precession in right ascension (degrees 

per day) 

Aa = correction in right ascension at time t 

for aberration (as defined on p. 47, 
"Supplement") 

In order to refer all data to the common reference time 
t^, it is convenient to rewrite (8) as follows: 
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(9) 



a(t) = a’Ct^) + a(t - t^) - a(t^ - t^) + Aa(t) 

= a'(t^) + a(t - t^) + AaCt) (10) 

where 

a'(t^) = a’(t^) - a(tj^ - t^) (11) 

Inserting (10) into (6) , 

GHA = GHA^^ + a(t - t^) - a’ (t^) - a(t - t^) - Aa(t) (12) 

= [GHA^^ - a’(t^)] + (a - a)(t - t^) - Aa(t) (13) 

Or, 

GHA = a^g + a^g(t - t^) - Aa(t) (14) 

where 

=*18 = GHA^o - 

- GHA^^ - a(tj^ - t^j 

- a' (t^) + a(t^ - t^) (16) 

and 

a^g = a - a (17) 

Now Aa must be developed into a usable form. The equation 
on pages 47 and 48 in "Supplement" [Ref. 8] shows that aber- 
ration at time t is given by 

= -k sec 6' [cos a’ cos X cos e + sin a* sin X] (18) 

where X is the longitude of the mean sun measured along the 
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ecliptic from y > e is the obliquity (angle between the mean 
planes of the ecliptic and equator) , k is the constant of 
aberration defined on page 48 in "Supplement , having a 
value of 20.47", and a' and 6' are the true right ascension 
and declination, all values being at time t. Selecting t^ 
as a convenient time for determining the mean place of the 
star and the obliquity, and regarding a', 6' and e as 
constants over the interval from t^^ to t, equation (18) 
becomes : 

Aa = - S6C cos cos e^^) cos X 

- sec sin ap sin X (19) 



= a^gcos X + a^^^sin X 



where 



a 



10 



20.47 

3600 



sec 



cos cos 



^1 



a 



11 



20.47 ^ 
3600 



4 



sin 



(20 



( 21 ) 

( 22 ) 



The factor 3600 enters because Aa will be computed as 
fractions of a degree whereas k was giv'^en in arc seconds. 
Equation (20) can be rewritten as 

Aa = a ^2 sin (X + a^^) (23) 
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10 



+ a 



11 



(24) 



and 



13 



tan’^ Ci^) 
^11 



(25) 



The longitude of the mean sun increases linearly with 
time according to the equation 



X(t) = X(t - t^) + X^ 



(26) 



where X^ is the known longitude at some convenient reference 
time tj^. As has previously been done, it is best to refer 
to a common reference time t^ by rewriting (26) as follows; 



X(t) = X(t - t^) - X(tj^ - t^) + X^ 



= X(t - t^) - X^ 



(27) 

(28) 



where X is defined as follows: 
o 



^o ^1 ’ ^^^1 ’ ^ o ^ 



(29) 



Substituting (28) into (23), the results are 



Aa = a, ^ sin [Xft - t ) + X + a^„] 
12 ^ 0 0 13 -^ 



a^2 sin [X(t - t^) + a^^] 



14- 



(30) 

(31) 
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where is defined as follov/s : 



a, . = X + ai , 
14 o 13 



( 32 ) 



Using (31) in (14) the substitution yields 
GHA = a^g + a^g(t - t^) - a^^ ^14^ 

Because a^^g includes earth rate, roughly 360° per day, and 
that t - t^ may be on the order of 30,000 days, to preserve 
accuracy in the computation it is v^rorthwhile to separate 
a^^g into two parts, one of which, ^ 21 ’ 360° and the 

other ^ 22 ’ remainder. Let 



^19 ^21 ^22 



(34) 



where 



= 360 (35) 

and 

^22 ^ ^19 ” 

In similar fashion, denote the number of days from the 

reference t to the time t to be D.d where D is a whole 
o 

number of days since 0000 Nov 29, 1971, and d a fractional 
part. That is 



t-t^=D.d (37) 

o 

Then it is obvious that the product of the angular rate times 
the time yields an angle modulo 360° which can be expressed 
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as follows : 



(360 + a22) (D + d) 


(38) 


360D + a22(D + d) + 360d 


(39) 


a22(D + d) + 360d 


(40) 


a22(D.d) + 360d 


(41) 



Substituting this result and (37) into (33) yields 



GHA = a^g + + 360d - ^14^ (^^2) 

The constants a and 6, precession in right ascension and 
declination, for a particular star are obtained using the 
formulas on page 39 in "Supplement," except that the values 
of the constants have been taken for epoch 1975 from the 
American Ephemeris and Nautical Almanac [Ref. 10] . Conver- 
sion to degrees per day is 



a 

6 



3'600 X 365.24219418 [46.1060 + 20.0404 sin a| tan <S|] 
TbOO X 365.24219418 [20-0404 cos a^] 



(43) 

(44) 



This now completes the derivation for the Greenwich 
Hour Angle of a star at time t. 

Turning attention to the derivation of declination, it 
is noted that the derivation is similar to those just 
carried out. Declination of a star changes because of 
precession and parallax. Letting 6 be the apparent 
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declination at time t, 6' be the true declination at some 
specified time, 6 be the rate of precession of this star 
in declination, A6 be the correction to declination for 
aberration, the result is 

6(t) = 6'(t^) + 6(t - t^) + A6(t) C45) 

where t^^ is the reference time of which the true right 

ascension and declination are known. The numerical value 

of 6 is given in (44) . This form is converted to the 

standard reference time t by the same method which has 

o 

been used above. The result is 

«5(t) = 6'(t^) + 6(t - t^) + A6(t) (46) 

where 

«S'(t) = 6'(t3L) - 6(t^ - t^) (47) 

From pages 47 and 48 in "Supplement" the time behavior of 
the correction for aberration is 

A6 = -k(sin cos - cos sin a| sin 6^) cos A 

-k cos sin sin A (48) 



— 


a^^^ cos A + sin A 


(49) 


where 








^15 


20.47 

3600 


(sin e| cos - cos a| sin aj sin 6|) 


(50) 


^16 


20.47 

3600 


cos sin 


(51) 
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The constants and a^^ are in degrees. Simplifying further 



- 1 i; 

A6 = a, - sin (X + tan — 

^16 



(52) 



where 



‘17 



2 ^ 2 
^16 ^15 



(53) 



Substituting (52) and (28) into (46) , 



6(t) = 6'(t^) + 6(D.d) + a^y sin[X(t - t^) + a^^^] (54) 



o" “20- 



where 



- 1 ^15 

a„rt = X + tan 

20 o a^^ 



(55) 



and 



t - = D.d 



(55) 



The final equations. of importance which lead to GHA and 
declination are (42) and (54) . Note that eight star- 
particular constants are involved in these two equations, 
plus X, a constant which appears in all calculations of 
GHA and declination of the stars. 

The equations involving reference times have been set 
up in such a way as to maximize freedom in their selection. 
Reference time t was selected to be 0000 Nov 29, 1971, to 
facilitate the calculation of time in days with a minimum 
number of entries as suggested by Ref. 8. 
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The reference time is any convenient time, presumably 
in the recent past, at which the GHA of Aries is precisely 
given in some reference tabulation. The time t^ is referenced 
to true right ascension and declination as listed in some 
reference almanac. The true right ascension and declination 
at a'(tj^) and 6 ' (t^) can be obtained from the apparent 
values a(t^) and <S(t^) by setting t = t^ in equations (8) 
and (45) and solving for a' and 6'. 

The reference time t^ is also the time, again in the 
recent past, at which obliquity e^, aj*, and are 
accurately known. The reference time for X^ is such that 
the accuracy of the sun's mean longitude is known from some 
almanac data. It will minimize the errors due to (a) eccen- 
tricity of the earth's orbit, and (b) the leap year cycle, 
if t^ is chosen at or near the time of earth perihelion, 
in the midpoint of a leap year cycle. 

It should be noted that t^ for the known value of GHA 

of Aries may be different from the t^^ for the known values 

of true right ascension and declination. Likewise t^ may 

be yet another time for the known sun's mean longitude. 

Since all the values (constants) associated with t^ are 

retarded to t„ no difference in t, is objectionable, 
o 1 

B. IMPLEMENTATION 

Tlie program presented in Appendix A is the result of 
several conversions. The first conversion consisted of 
taking algorithms expressed in Ref. 8 with the HP- 65 
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programs to produce a debugged FORTRAN IV program. The 
HP-65 programs deal with degrees as a unit of angular measure 
vice radians. A conversion of all data from degrees to 
radians was made for use in the FORTRAN IV and BASIC pro- 
grams. General routines 0010 and 0020 provide for I/O 
conversion of radians to degrees, and hours, minutes, and 
seconds to hours (and vice versa) . These conversions were 
needed for ease of data entry and data output. 

The HP- 65 programs are compact due to the limited 
number of registers used for intermediate storage. This 
efficiency of memory was carried through to the FORTRAN IV 
programs and then to the BASIC programs for the routines 
that determine the GHA and declination of the star and the 
subroutine that determines the azimuth (Zn) and distance 
(subroutine 0300 and 0400, respectively). It was felt that 
this level of efficiency need not be continued in the 
remainder of the program due to the need for self -documentation . 

FORTRAN IV was selected as the intermediate step in the 
conversion due to the amount of core memory allocated to 
the FORTRAN compiler. The BASIC interpreter is limited in 
memory on the IBM 360 installation and would not hold the 
completed BASIC program. 

The reduced precision of the FORTRAN IV version of the 
program (as it was not programmed in double precision) was 
attributable to the internal formatting of the variables. 

High precision was not of great concern at this stage of 
development . 
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As the goal of this process was to provide a program 
that would implement the algorithms necessary to complete 
a celestial fix on a microprocessor, the FORTRAN IV version 
was converted to BASIC. The BASIC program allows for much 
greater precision in computation as the values are expressed 
in 32 bit floating point, when interpreted by the available 
BASIC interpreter in the Microcomputer Laboratory at the 
Naval Postgraduate School. The BASIC interpreter at the 
W. C. Church Center was not of sufficient size to run the 
entire program. Therefore, the main subroutines and all 
required support routines were run and debugged separately, 
and required information was transferred manually from rou- 
tine to routine. Because the BASIC language does not pro- 
vide for scopal variables and due to the limited types of 
identifiers available much of the self -documentation of the 
FORTRAN IV program was lost. Table I correlates variables 
in subroutine 0100 in the program listing (Appendix A) to 
those variables previously discussed during the theory 
presentation . 

Subroutine 0100 computes the eight constants associated 
with a star and places them in array C (elements C(2)-C(9)). 
Element C(l) is the rate of increase of the mean sun longi- 
tude (common to all stars) . 

In order to organize the variables as much as possible 
a scheme was invoked as explained below. Since there were 
a number of subroutines in the program, a method of passing 
variables was constructed. The variables R1-R9 and P1-P9 
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TABLE I 



Identifier Correlation of Variables in Subroutine 
0100 and Equation Variables 

Subroutine Variable Equation Variable 






C(l) 


X 


C(2) 


^20 


C(3) 


^17 

e 


C(4) 


6 


C(5) 


6'( 


C(6) 


^22 


C(7) 


^12 


C(8) 


00 


C(9) 


^14 


13 


ai 


14 


O 


IS 


a 


16 


&19 


17 


^10 


19 


^11 


JO 


^15 


J1 


^16 


C5 


^1 
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are passing variables. Variables used in three or more 
non-general routines were considered as global and were 
selected to be as self-documenting as possible. For example, 
the two constants for each star are placed in the array 
S(X,X). The constants used by routine 0300 to determine 
GHA and declination which are created from the S array are 
placed in the C vector. All other constants were repre- 
sented by the format C1-C9. 

Several routines use intermediate variables while calcu- 
lating the object variables. These variables are repre- 
sented by 11-19, J1-J9 and M1-M9. All logical values are 
represented by L1-L9. Note that this scheme was invoked 
to eliminate problems created by not having scopal varia- 
bles. This must be considered in the selection of a 
language for other peripherals used in the system. Points 
to be considered in selecting a final programming language 
for such a system and the peripherals attached to it are 
discussed in Section VII. 
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VI. THE INFORMATION DISTRIBUTION SYSTEM 



A. OBJECTIVES 

This subject will primarily be concerned with the 
integration of a general purpose computer into the existing 
shipboard information flow. The "system" in the spotlight 
is the navigation suite aboard U.S. Navy ships. The phil- 
osophy for this approach is simple and is best stated by 
Joseph Weir [Ref. 1] : 

"Good systems can be designed only in the presence of 
potential users and applications. By using the system 
in numerous applications and assiduously watching its 
performance, it gradually can be adopted to the (actual) 
needs of the application and the user. Ingenious ivory 
tower solutions frequently do not v;ork or they solve 
problems that don't exist." 

A close examination of this statement while relating 
it to the navigation suites \\^ill show that in fact the 
"system" has been designed both as far as it need and can 
be with the existing equipment. Note that the "system" 
comprises geographically dispersed equipment and telephones 
for communication links. 

The application and needs of the navigation system have 
been defined as will be discussed later in this section. If 
all this is true then why rock the boat by imposing a 
system, apparently complex, on an existing "system"? 

Discussion with crewmen aboard the USS Hancock indi- 
cated that the existing system is cumbersome, subject to 
human frailties and is not timely enough for the devices 
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on board. OMEGA radio navigation sets are capable of 
providing information for navigational fixes every few 
seconds, yet this information is only used hourly. In 
this case the communication link is the "weak” link, and 
to eliminate this weak link an information distribution 
system must be imposed. The distribution system must 
appear translucent and straightforward to the user, for as 
Weir mentions, 

"The more complex the system, the more likely it is 
that it won't get used or will be misunderstood or 
misused." [Ref. 1] 

All systems have sensors and effectors, and some 
communication constraints placed on them. Although most 
of these sensors for the navigation system in fact exist, 
there are always new devices being created. 

For example, ten years ago operational OMEGA navigation 
sets did not exist. If a system had been created and imple- 
mented on the ships at that time without forethought with 
regard to new devices, then the information system would 
not have been able to integrate the new device into the 
network . 

In order to fully utilize an information distribution 
system there must be "terminals" with which to input informa- 
tion and query the system. For example, suppose the bridge 
wants to know the average RPM of the four propulsion shafts. 
Engineering would have a terminal (either automated or manual) 
that would accept certain required data, e.g., tachometer 
inputs, transform it into its varied forms and relay the 
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data to the central processor. The processor would then 
transfer the required information to the terminal at the 
bridge . 

B. DESIGN PHILOSOPHY 

The question might be asked, "Why try to integrate an 
information distribution system with the existing periph- 
erals?" The answer should be obvious - -cost of implementa- 
tion. Existing equipment must be fully exploited until 
it is not wise to continue using it. 

By integrating existing equipment into a network instead 
of creating all of the component parts of the system anew, 
the cost of research and development of all the new sensors 
and effectors in the system can be elided [Ref. 2] . Of 
course, the use of current equipment will demand the design 
of interfacing equipment such as analog to digital converters, 
digital to analog converters, fault sensing networks, etc. 
However, the cost of these units should be relatively inex- 
pensive when compared to the redesign of any one sensor or 
effector in the system. The application of the same item 
(A/D or D/A, etc.) to many peripherals is a cost savings 
when design efforts are considered. 

As distribution of existing data is the concern of the 
central system and since there is a multiplicity of like 
data, the question of which source's information should be 
distributed must be handled. A rather trite answer is 
"the best." Again a question, "Which sensor is providing 
'the best' information?" 
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A definition of the term "best” is recursive and must 
include words like timeliness, accuracy and statistical 
validity. As each sensor in the system is operating at 
pre-system design rates, some periodic (auto track LORAN C) 
and others at non-periodic rates (celestial sightings) , 
timeliness takes on a meaning associated with the quality 
of the source. For example, if the only peripheral capable 
of input to the system is the celestial navigation unit, 
an output of ship's position would only be considered timely 
if the last celestial position had been updated with inputs 
from the "best" heading and speed information to provide an 
estimate of the current position. 

The word "accuracy" when associated with "best" means 
most accurate on a timely basis. A position based on a 
ten-hour old OMEGA position is usually not as accurate as 
a five-minute old celestial position. 

This is all well and good, yet how are these descriptors 
compiled? A central machine is obviously the device that 
will have to provide the software to meet problems based 
on statistical properties. 

Statistical properties of each measuring device must 
be taken into account. The proper error probabilities must 
be applied to each input; then based on the result of the 
combined statistics, the "best" source is selected to repre- 
sent that class of input. It is this source's information 
(corrected for time) that must be distributed. 
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Along with statistical constraints imposed on the data, 
a status check must also be im.posed on the device. If 
the data is apparently correct yet the device has failed, 
then the data is obviously of no use. The status determina- 
tion system must be as nearly fail-safe as possible. Since 
this paper is proposing a centralized star network vice a 
fully connected or ring network [Ref. 2] , all the status 
checks must be initiated at the central computer (see 
Figure 6). In order to design a fail'safe status monitor, 
the following considerations must be employed: 

1. Notice of failure of the central computer must be 
distributed to all peripheral stations. 

■ 2. Failure of each peripheral must affect only that 

peripheral and, at worst, place the system in a degraded 
mode . 

3. Status of selected peripherals might have to be 
recorded for legal records. 

4. Failure of the peripheral must not cause failure 
of the status check mechanism. 

There are several reasons for designing an information 
distribution system. Many operations aboard a ship are 
dependent on the same source of information. For example, 
the Combat Information Center and radar stations need true 
course and speed information. These are parallel operations 
and their information must be based on the same source. A 
distribution system would be the answer to this problem. 
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All of the discussion thus far has been about distribution 
of information gathered from sensors, so where is all this 
fine information used? An example may best describe the 
answer. Suppose the navigator updates the current ship’s 
position with the celestial navigation computer. If the 
Captain desires to proceed to a desired point, the point 
could be placed into the system and the system would pro- 
vide proper steering information (great circle steering) 
to the Helmsman, or this data could be transmitted to an 
auto-pilot upon request. 
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VII. LANGUAGE SELECTION 



The subject o£ language selection is placed at this 
point in the paper due to its applicability to both the 
celestial navigation microcomputer and an information 
distribution system as described within. The considerations 
in selecting a programming language for these two devices 
have much in common, the most important of which is the 
microprocessor . 

There are two main considerations in selecting a pro- 
gramming language for a microcomputer. The first and most 
important consideration is that of efficiency; efficiency of 
both memory utilization and program execution speed. The 
second point concerns the cosmetics of the language (i.e., 
is the language easy to read and is it self -document ing?) . 

In consideration of the first property, that of 
efficiency, the language must be chosen in the light of 
cost/availability of volatile memory. Even though the 
addressing schemes of some microprocessors provide for 
large memories there may be physical and/or cost limita- 
tions placed on the microcomputer design that make very 
efficient use of the available memory imperative. For 
peripherals operating with unmodifiable code the number of 
memory locations can be optimized and placed with a high 
degree of efficiency. For processors using programs v/hich 
have subroutines witli arrays of data created during the 
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execution o£ the program, the placing o£ these arrays 
on a metal mask chip or £ixing and reserving the addresses 
o£ these arrays is an ine££icient use o£ memory. 

The best v^^ay to program in an e££icient manner is to 
use block structured programming techniques utilizing 
"scopal" variables. Block structured languages like Algol 
60 and PLM [Re£. 12] allow £or scopal variables to include 
arrays which provide £or e££icient use o£ random access 
memory and reduce the amount o£ programmable/read only 
memory. By creating arrays as needed, common random 
access memory may be used very e££iciently allowing £or 
lower investment in space and cost. 

The second area o£ concern, that o£ cosmetics, is a 
growing concern with all languages. It may be a seemingly 
unimportant aspect until the prospects o£ non-volatile 
runtime modi£iable LSI memories (e.g., £erroelectric 
memories [Re£. 15]), multiplicity o£ similarly constructed 
peripherals and common maintenance support are considered. 
The celestial navigation microcomputer program bases 
all o£ the calculations £or geographical positioning 
o£ the body constants which vary £rom year to year. 

In order £or the device to enjoy extended longevity, 
the constants must be alterable. The language must 
allow £or the documentation to help provide this support. 
An assembler program is not a sel£- documenting 
program no matter how e££icient it is. Its logic and 
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meaning are often forgotten with the loss of the programmer. 
Indeed, the programmer may lose the logic to his own program 
if he has been away from it for a length of time. For 
example, the statement below is obvious in its intent: 

NO$YEARS = NO$DAYS / 365. 

Yet even though the next statement provides for the same 
efficiency of calculation it does not connote as clear a 
meaning : 



LIT DAYS; LIT 365; DIV; LIT YRS ; STO; 

The assembler language form of this statement loses the 
easily conceived meaning of the first statement. 

Though the celestial fix program in Appendix A was 
written in the BASIC language due to the precision attainable 
through its use, the program itself would not be under- 
standable \^rithout the comment statements and extra 
documentation. 

Specifically, there are primarily two high level 
languages available to the microprocessor programmer; 

BASIC as implemented at the Naval Postgraduate School and 
PLM. BASIC, as the name implies, provides for the basic 
needs of a programmer. BASIC has subroutines, iterative 
and logical capability, the full range of mathematical 
manipulators and limited I/O. A quick glance at the program 
in Appendix A provides the evidence of a major drav\fback: 
lack of variable scope and identifiers. There are no local 
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variables; all variables are global. This property coupled 
with the limited identifier schemes causes a prograramer of 
a large program to run out of variables quickly. His atten- 
tion is turned to efficient use of identifiers instead 
of efficient programming techniques. 

PLM, on the other hand, provides for scopal variables 
and identifiers that allow for self-documentation . The 
following statements are equivalent; however, the second 
connotes more meaning: 

LET J1 = G(Q) - G(R) 

DIFF = GHAl - GHA2. 

The array used in the determination of the three 
intersections used to calculate the fix (subroutine 0600 . . . 
ARRAY V(6,6)) is a temporary array. In BASIC the array is 
permanently allocated space. Using PLM and scopal variables 
the array would be created by the program during execution 
(indeed, could be of variable size) and be destroyed or 
written over upon exit from the subroutine. There are a 
few drawbacks with the use of PLM. The primary drawback 
in the execution of the program in Appendix A is the word 
size. Sixteen bits do not provide the accuracy needed for 
the celestial navigation problem. Though this problem can 
be solved, it is not easily solved. 

Special mathematical routines like trigonometric func- 
tions are not available in PLM, yet they can be programmed. 
The point is that PLM can be tailored to the needs of the 
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programmer with a little extra work while the BASIC language 
cannot be tailored to the user at all. PLM provides a 
much more versatile language yet the detail required in 
programming is greater than that of BASIC. 

A final comment on the selection of a language deals 
with the multiplicity of processors on an information 
distribution system. If each processor were programmed at 
a different level with a different language the maintenance 
and design problems would be horrible. For example, suppose 
one processor is to provide a second processor information. 
If the first processor is program.med in BASIC then it would 
probably transmit its data in 32 bit floating point form. 

The second processor (programmed in PLM) would have to have 
an I/O routine to convert that processed data to a local 
format. If there were a third language then the PLM- 
programmed processor may need yet another I/O processing 
routine. This situation is the grounds for selection of a 
single programming language for a system of microprocessors. 
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VIII. CONCLUSIONS 



Currently, the proposed computer is capable of developing 
a navigational fix (without use of tlie Nautical Almanac and 
H.O. Publications) utilizing only the 57 navigational stars. 
Because of this limitation it is obvious that numerous and 
known applications are yet to be incorporated. Among these 
include methods which employ the attributes of the sun, moon 
and planets so that positioning with these bodies can be 
performed. Considering functional aspects for future 
development of programs, suggested programs are closest 
point of approach (CPA), great circle calculations, star 
location utilizing estimated altitude and azimuth and fix 
triangulation routines. 

The sun, moon and planets present special problems 
because of their proximity to earth. They do not appear 
fixed in the sky as do stars. Stars do mot^e across the 
sky due to the earth’s rotation but their great distances 
impart a constant relationship to their neighboring stars. 

The planets and moon, being closer to earth, possess no 
fixed relation in the evening sky with the stars. 

The size of the sun and moon requires a special con- 
sideration termed semi-diameter. More accurate sight reduc- 
tions can be performed by aiming the sextant at the upper 
and lower limb of these bodies; corrections for the radius 
of the body can be programmed. 
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For lunar observations horizontal parallax is a 
parameter required to be considered. It is the apparent 
difference in the position of a body when viewed from two 
different locations. Horizontal parallax is significant 
for the moon because it is the closest body to earth. 

The celestial lines of position are included to give 
the navigator the option of plotting LOP's. This method 
should provide the navigator with an azimuth, a distance 
and a toward/away indication in reference to the assumed 
position of the observer with which the LOP should be 
plotted. This is the backup routine to the fix calculation 
algorithm. 

The proposed celestial navigation computer was used as 
a vehicle for the approach to a design of a dedicated micro- 
computer. Other devices of similar construction used as 
peripherals in the information distribution system would be 
approached in the same manner. Approaching the design of 
each device on its own merits will allow for the system to 
grow in a manner that "fits" the situation, i.e., a single 
"ivory tower" system would not result. 

Indeed, the development of the central processor 
(assuming enough ports are planned for) would grow with the 
system as each new unit would simply add length to the polling 
routines and interrupt structures. Following along with 
Weir’s natural growth concept, each processor's programs 
would be designed and built (or at least the algorithms be 
provided for) as needed by cognizant individuals, from 
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current literature, as done in the case of the celestial 
navigation computer. The only restraint placed on the 
design/designer would be that of language and interface 
specifications. As delineated earlier, the authors believe 
that of the two languages examined PLM would be the most 
efficient language and would provide for much needed 
self -documentation . 

The central processor in the information distribution 
system may very well be a minicomputer or a microcomputer 
depending on its functions. It need only be a distribution 
device designed along the same lines as the peripherals 
it supports. The level of computation provided by the 
central processor must be established early in the design 
phase because the peripherals would need to "take up the 
slack" for the computational capabilities not provided by 
the central processor. For example, the averaging of the 
revolutions per minute may be accomplished at the peripheral 
(which is recommended) or at the central processor on the 
request of the inquirer. Again, keeping this device 
simple (only the necessary switching and monitoring routines) 
will provide for more independence and faster information 
transfer . 
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31 M B(2),L(6),E(S),V(6,6),D(3)fH(3),G(5) 

DIM C(9),Z(3)tA(3),0(3)tS(2,57) 

REM DEFINE LOGICAL VALJES AND CONSTANTS 



APPENDIX A 



CELESTIAL FIX PROGRAM LISTING 
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APPENDIX B 



SAMPLE PROGRAM RUN 



ENTER DAY 
20 

ENTER MONTH 
6 

ENTER YEAR 
75 

ENTE R ALTI T JOE { FT MSL ) 

60 

ENTER COURSE (D.MS) 

270 

ENTER SPEED (KNOTS TRUE) 

5 

ENTER ASSUMED LAT (N/S) +/-(□. MS) 

40 

ENTER ASSUMED LONG(E/W) +/-(D.MS) 
15.2224 

ENTER NUMBER OF BODIES THIS FIX 
3 

ENTER STAR/BDDY NUMBER 
1 

ENTER GMT 3- DBS (H.MS) 

.56 

GMT (HRS)= .944443 

ENTER HS (D.MS) 

37.04 

ENTER IC (D.MS) 

0 

HA= 36.5709 

GHA OF BODY= 230.043' 

DEC OF BODY (D.MMSS) 

IN OF BODY (D.MMSS) 

DISTANCE TOWARD BODY (NM) 

BASED ON H: (D.MMSS) 



28.5728 

8D.423). 

8.25D71 

36.4854 



ENTER STAR/30DY NUMBER 
46 

ENTER GMT OF OBS (H.MS) 

1 

GMT (HRS)= 1 

ENTER HS (D.MS) 

48.39 

ENTER IC (D.MS) 

0 

HA= 48.3209 

GHA OF BODY= 19.1425 

DEC OF BODY (D.MMSS) 

ZN OF BODY (D.MMSS) 
DISTANCE TOWARD BODY (NM) 
BASED ON HC (D.MMSS) 



12.3436 

237.462 

-30.9003 

49.3303 



ENTER STAR/BODY NUMBER 
40 

ENTER GMT 0= OBS (H.MS) 
1.04 

GMT (HRS) 1.07778 

ENTER HS (D.MS) 

42.11 

ENTER IC (D.MS) 

0 

HA= 42.0409 

GHA OF BODY= 61 . 1234 
DEC OF BODY ( D.MMSS) 

ZN OF BODY ( D .MMSS ) 
DISTANCE tOWARD BODY (MM) 
BASED ON H: (D.MMSS) 



74. 1 533 
339. 153 
14. 4483 
41.4942 
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tNTER SHDT DF DESIRED FIX TIME 
2 

CORRECTED DISTANCE TOWARD BODY 
7.92129 

BASED ON CORRECTED HO (D.MMSS) 
36. 483A 

CORRECTED DISTANCE TOWARD 300Y 

-30. 9003 

BASED ON CORRECTED HO (D.MMSS) 
49.0303 

CORRECTED DISTANCE TOWARD BODY 
14.3302 

BASED ON CORRECTED HO (D.MMSS) 
42.0402 

FIX (LAT/LONG D.MMSS)= 

40.3313 15.3354 

TIME OF FIX ( .DAY)= 

PROGRAM COMPLETE 



1 



2 



3 



.D41S566 
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